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ABSTRACT 

The association of long-duration gamma-ray bursts (LGRBs) with Type Ic supemovae presents a challenge to 
supernova explosion models. In the collapsar model for LGRBs, gamma rays are produced in an ultrarelativistic 
jet launching from the magnetosphere of the black hole that forms in the aftermath of the collapse of a rotating 
progenitor star. The jet is collimated along the star's rotation axis, but the concomitant luminous supernova 
should be relatively — though certainly not entirely — spherical, and should synthesize a substantial mass of 
56 Ni. Our goal is to provide a qualitative assessment of the possibility that accretion of the progenitor envelope 
onto the black hole, which powers the LGRB, could also deposit sufficient energy and nickel mass in the 
envelope to produce a luminous supernova. For this, the energy dissipated near the black hole during accretion 
must be transported outward, where it can drive a supernova-like Shockwave. Here we suggest that the energy is 
transported by convection and develop an analytical toy model, relying on global mass and energy conservation, 
for the dynamics of stellar collapse. The model suggests that a ~ 10, 000km s _1 shock can be driven into 
the envelope and that ~ 10 51 erg explosions are possible. The efficiency with which the accretion energy is 
being transferred to the envelope is governed by the competition of advection and convection at distances ~ 
100- 1 , 000km from the black hole and is sensitive to the values of the convective mixing length, the magnitude 
of the effective viscous stress, and the specific angular momentum of the infalling envelope. Substantial masses 
of 56 Ni may be synthesized in the convective accretion flow over the course of tens of seconds from the initial 
circularization of the infalling envelope around the black hole. The synthesized nickel is convectively mixed 
with a much larger mass of unbumed ejecta. 

Subject headings: accretion, accretion disks — black hole physics — gamma rays: bursts — nuclear reactions, 
nucleosynthesis, abundances — supernovae: general 



1. INTRODUCTION 

A growing number of long-duration gamma-ray bursts 
(LGRBs) are being discovered in association with Type Ic su 
pernovae (|Galama et al.| [J998, 2000; Reichart 1999; Bloom 
Valle et al1|2 003 2006; Garnavich eTal 
2003 ; Hjorth et al. 2003 ; Kawabata et al. 2003 ; Stanek et al. 
2003 ; Matheson et al. 2003 ; Malesani et al. 2004; Campana et 
al. 2006; Mirabal et al. 2006; Modjaz et al. 2006; Pian et al. 
2006; Chornock et al. 2010; Cobb et al. 2010; Starling et al. 



the stellar envelope for an explosion, and that 56 Ni is synthe 
sized in th e wind to later produce an opt i cally-bright super 
nova (e.g., IMacFadyen & Woosley||1999| |MacFadyen||2003 
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r-\ Pruet et al. 



2003 2004; Kohrietal. 



2005[ ). The mechanics of 



201 1 ), yet on the basis of the non-detection of late-time radio 
emission in a sample of Type Ic supernovae, Podsiadlowski 



et 
than 



d] ( 12004} and |Soderberg et al.| ( |2006) inferred that less 
irilO% or all Type Ic supernovae are associated with stan- 
dard LGRBs. The process producing LGRBs an d their con- 
comitant supernovae remains a subject of debate (Woosley~&| 
Bloom 2006 and references therein). In the collapsar model 
for LGRBs (Woosley 1993), the gamma rays are produced in 
an ultrarelativistic jet launching from the magnetosphere of 
the black hole that forms in the aftermath of the collapse of a 
rotating progenitor. The jet is powered by a continuous infall 
and disklike accretion of the progenitor star's interior. While 
the collapsar model seems to successfully explain the power 
and duration of LGRBs, it is not clear at present whether it 
naturally gives rise to a supernova-like stellar explosion. It 
has been argued that a "wind" outflowing from the nonradia- 
tive parts of the collapsar disk may convey sufficient energy to 

1 Texas Cosmology Center, University of Texas, 1 University Station 
C1400, Austin, TX 78712. 

- NSF Graduate Research Fellow. 

3 Currently at Department of Astronomy, University of Toronto, 50 St. 
George St., Toronto, Ontario M5S 3H4, Canada. 



energy transfer from the disk to the supernova ejecta and its 
implications for nickel synthesis remain open problems and 
are the subject of the present study. 

Here we provide a crude assessment of the possibility that 
the accretion onto the black hole that powers the LGRB might 
also deposit sufficient energy in the progenitor envelope to 
produce a supernova. For this, the accretion energy dissipated 
near the black hole must be transported to exterior mass co- 
ordinates of the star. We hypothesize that the energy is trans- 
ported by convection to energiz e the ou tward moving shock, 
as was originally suggested by Narayan et al. ( 200 1 1> . Our 
two-dimens ional simulations of collapsar accretion ( |Lindner| 
|et al.|2010| , in which we simulated only relatively large radii 
(r > 500 km) and did not incorporate neutrino and nuclear 
physics, corroborate the crucial role of convection. Our ap- 
proach in the present work differs from existing assessments 
of the viability of collapsar supernovae (e.g., Kohri et al. 
|2005| l in that we attempt to sketch out the global structure 
of the flow by incorporating disk accretion, convective energy 
transport, shock dynamics, and long-te rm stellar infall (the 
latter having been studied in Kumar et al. 2008a b) in a sin- 

a I . ll_ __1_L 

gle toy model. Particularly important are the non-Keplerian 
nature of the accretion flow and the precise form of the vis- 
cous torque in a pressure-supported regime; the latter aspects 
seem to have been neglected in existing treatments of collap- 
sar disks bu t came to light in the 2.5D simulations of Lindner 
|eTaT1 ( |20T0T ). 
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The toy model allows us to investigate how the shock ex- 
pansion interferes with the rate with which the infalling ma- 
terial is accreting onto the black hole. This differs from 
the approaches that implicitly postulate, as in the advection- 
dominated accretion flow (ADAF) paradigm (e.g., Narayan & 
|Yi|1994|[l995j|Blandford"&Be gelman 1999), that the inflow 
occurs in the equatorial region and is spatially separated from 
a non-interfering wind that carries mass and energy into an 
axial outflow cone. Our treatment is similar in spirit to the toy 
model constructed by Janka (2001), following in the footsteps 
of Bethe (1990 1993a b 1995 1996 1997), to assess con- 
ditions for shock revival by neutrino heating in core-collapse 
supernovae; the key differences include the central roles of ro- 
tation and convection in collapsars, and comparatively lower 
mass accretion rates, longer time scales, and lower fluid den- 
sities in the latter systems. 

This work is organized as follows. In Section[2]we discuss 
the formation of an accretion shock and the post-shock con- 
ditions immediately following shock formation. In Section [3] 
we analyze the structure of the inner accretion flow and esti- 
mate the luminosity that convection can transport toward the 
shock wave. In Sectionfflwe impose global mass and energy 
conservation and in Section [5] we utilize these to so estimate 
the shock expansion velocity and the total energy deposited 
in the stellar envelope. In Section[6]we discuss the prospects 
for 56 Ni production in our model. In SectionMwe present our 
conclusions and briefly discuss some implications. 

2. THE ACCRETION SHOCK 

The formation of the black hole in stellar core collapse 
might be preceded by a failed explosion resulting from a 
bounce and subsequen t shock revival by neutrin o heating 
( |Bethe & WiIsonj[T9"85| see also |MacFadyen et al.pOOl) , or 
by an ejection of a magnetic field from a magnetized proto- 
neutron star (e. g., Bisnovatyi-Kogan 1971, Wheeler et al. 



2000; Thompson et al. 2004; Bucciantini et al. 2007, Burrows 
et al. 2007 ; Dessart et al. 2008 ). The magnetic outflow may be 



too axially collimated to produce a standard supernova explo 
sion (Bucciantini et al. 2008 2009). Here we ignore the pos- 
sibility of any type of explosion preceding the collapse into a 
black hole, and assume that during the first few seconds from 
the collapse of the stellar core, an unshocked stellar envelope 
accretes quasi-radially onto the black hole. 

For a few solar mass black hole, direct unshocked accretion 
through the event horizon of the black hole is possible dur- 
ing an initial interval measuring in the tens of seconds, while 
the specific angular momentum of the infalling shells is in the 
range I < 10 16 (M BH /M Q )cm 2 s" 1 (M BH is the black hole mass 
and the critical specific angular momentum decreases with in- 
creasing black hole spin for prograde accretion). Because the 
specific angular momentum of the innermost stellar shells in- 
creases outward, unshocked accretion takes place as shells at 
mass coordinates ~ (3 - 6)M Q arrive near the black hole in 
the unmixed and fully mixed pre-supernova models (see, e.g., 
Heger et al. 2000 2005; MacFady en et al.|[200T]|Petrovic et 



al. 2005; Woosley & Heger 2006). When the specific angu- 



lar momentum arriving near the innermost stable orbit around 
the black hole becomes comparable to the angular momentum 
of the innermost stable circular orbit (ISCO), the accretion 
occurs through a "dwarf" or "mini"-disk dLee & Ramirez- 
Ruiz 2006; Zalamea & Beloborodov 2009); during this transi- 
tional period, the black hole could acquire rapid rotation and 
accretion takes place at the rate M ~ (0.1 -0.2)M Q s _I . At 
some time f s h that depends on the stellar mass and rotational 



profile, the flow crossing the ISCO becomes subsonic and a 
quasispheri cal accretion shock wave fo rms aroun d the black 
hole (e.g., MacFadyen & Woosley 1999 ; Lee & Ramirez-Ruiz 
2006; Nagataki et al. 2007 ; Lopez-Camara et al. 2009; Lind- 
ner et al. 2010). 



The initial radius of the shock, r s h, is just larger than 
Isco ~ 5 -50km, where the latter depends on the mass and 
angular momentum of the black hole. If the shock proceeds 
to travel outward, the rate wit h which material ac cretes onto 
the black hole drops rapidly ( jLindner et al.|20K)) , e.g., by a 
factor of 10 or more during the first second from the inception 
of the shock. Dynamics of the shock wave is governed by the 
rate of stellar envelope infall and the conditions in the rotating 
downstream fluid. Besides the heating at the accretion shock, 
the fluid is heated by the dissipation of magnetohydrodynamic 
turbul ence driven by the m agnetorotational instability (MRI, 
e.g., [Thompson et aT. 2005, and references therein). The fluid 
cools by neutrino emission, primarily through pair annihila- 
tion and pair capture onto nucleons (the Urea process). The 
fluid also cools through the disintegration of nuclei into he- 
lium and free nucleons as it passes the shock and accretes 
toward the black hole. The latter process can be reversible, as 
the energy consumed in disintegration can be recovered if the 
free nucleons and helium end up getting transported by con- 
vection ("dredged up") to radii with lower entropies where 
they can recombine into heavier nuclei. 

The character of the flow is sensitive to the relative mag- 
nitude of the cooling and heating rates. When the cool- 
ing is comparable to the heating, the flow collapses into a 
rotationally-supported, neutrino-cooled accretion disk (e.g., 
MacFadyen & Woosley 1999; Popham et al. 1999; Narayan 
et al. 2001 ; Di Matteo et al. 2002; Kohri & Mineshige 2002 
Janiuk et al. 2004; Setiawan et al. 2004; Kohri et al. 2005 
Lee et al. 2005 ; Chen & Beloborod ov|20071|Kawanaka&Mi-| 
neshige 2007). The flow accretes onto the black hole through 



a thin disk when the accretion rate is higher than a minimum 
value that depends on the viscous stress-to-pressure ratio a 
and the black hole's spin parameter a; this is because disks 
with a larger a are less dense and cooler. In disks that are op- 
tically thin to neutrinos, for a given a, neutrino cooling dom- 
inates disk thermodynamics at relatively high accretion rates 
and small radii; at lower accretion rates or larger radii, the 
flow becomes g eometrically thick a nd non-radiative. For ex- 
ample, [Chen^&Beloborodov] ((2007J) found that for a = 0.01, 
a neutrino cooled disk can be present around a Mbh = 3M Q 
black hole for M > (10" 4 - 1O" 3 )M s _1 , but for a = 0.1, the 
accretion rate must be M > (0.02-0.1)M Q s" 1 for a thin disk 
to be present. ' 

The true value of a in the regime in which pressure and 
centrifugal forces are of the same order, and where the fluid, 
as we shall see, is convective, is not known. If the flow is con- 
vective, then the convection may contribute to the buildup of 
the magnetic stress ( Balbus & Hawley 2002 Igumenshchev 

3J), 



|2002| |Igumenshchev et al.|2003[ jChristodou ou et al.|200; 



1 If the volumetric neutrino cooling rate is Q v <x [/* T B , where (A,B) = 
(0,9) for pair annihilation and (A,B) = (1,6) for Urea, and if, for the pur- 
pose of illustration, radiation and low-density relativistic pairs dominate the 
pressure at the thin-to-thick transition radius r v and the gravitational field 
is Newtonian, then it is straightforward to show that for fixed M, we have 
r„ <x a 2 (4A+B)H32-l2A+5B) ^ ^j^ gives r„ <x a -18 / 13 for pair annihilation and 
r v <x aT 2 for Urea. The value of a must be smaller than a critical maxi- 
mum value if r„ is to be larger than nsco, as required for neutrino-cooled, 
think-disk accretion. 



and this may motivate a large value of a and a large criti- 
cal accretion rate required for the presence of a thin, neutrino 
cooled disk. If a ~ 0.1, which is the value we take as the 
fiducial for what follows, then the very initial accretion rate 
drop following shock expansion already brings the accretion 
rate below the critical rate for efficient neutrino cooling, and 
the flow is nonradiative and geometrically thick at all radii. 2 
The toy model that we present below will be restricted to this 
non-neutrino-cooled regime. 

It may be worth noting that the picture in which a super- 
sonically infalling flow passes an accretion shock, becomes 
predominantly rotationally supported, and proceeds to ac- 
crete onto a central compact object, either through a neutrino- 
cooled thin disk, or through a nonradiative thick disk, res- 
onates with the work examining the post-supernova fallback 
onto a neutron star, or examining the Bondi-Hoyle accre- 
tion onto a neutron star embedded in a common envelope 
(see, e.g., Chevalier 1996; Brown et al. 2000 and references 
therein). Because the characteristic accretion rates in these 
contexts, which are ~ lM^yr' 1 , are orders of magnitude be- 
low those anticipated in collapsars, the maximum values of 
a for which neutrino-cooled disk solutions exist are much 
smaller than those in collapsars. Furthermore, photon diffu- 
sion may be relev ant in the fa llback and common-envelope 
contexts (see, e.g., Blondin|198 6), but in collapsars, complete 
photon trapping is a safe assumption. 

We work under the assumption that the accretion shock re- 
mains quasi-spherical as it traverses the star, and thus, that the 
thermal "wind" produced in the inner accretion disk remains 
trapped within the surface of the shock so that the wind's en- 
ergy is distributed quasi-spherically beh ind the sho ck, as seen 
in idealized 2.5D simulations ( [Lindner et al.|2010) . We can 
allow for the possibility that a collimated electromagnetic out- 
flow distinct from the thermal wind, such a s a jet enveloped by 
a cocoon of shocked stellar fluid (see, e.g., Zhang et al. 2003 



2004 
eTaT 



Zhang & MacFadyen 2006 1 Morsony et al. 2007 ; Wang 
2008 ), is present along the axis of rotation; our analy 



sis should be construed as applying to the equatorial region 
not occupied by the jet. Regardless of the presence of the jet, 
on time scales much shorter than the free fall time from the 
surface of the star, the axial "funnel" region is not empty, and 
remains overpressured either by the freely falling axial low- 
angular-momentum material, or by the jet's hot cocoon. Ther- 
mal outflow from the predominantly rotationally-supported 
central accretion flow launches at oblique angles from the 
surface of the disk, following "gyrentropes" (the surfaces of 
approximately constant angular momentum, Bernoulli func- 
tion, and entropy; see Blandford & Begelman 2004), as is 
evident in numerous simulations of radiatively-inefficient ac- 
cretion flows in regio ns in which the magnetic field is not dy- 
namically important (Stone et al. 1999; Igumenshchev et al. 
2000 ; Abramowicz et al. 2002 ; Hawley & Balbus 2002 ; Proga 
& Begelman 2003; Prog a et al.||2003||Lindner et al -p-OWf 
The thermal wind then mixes with the post-shock fluid; in 
this sense, we think of the oblique thermal wind as a form of 
convection (or, more adequately, stochastic circulation) with 
an effective mixing length that can be large and need not be 



2 In our companion 1 .5D numerical simulations {Lindner et al.|201 1| , we 
estimated the vertical-pressure-scale-height-to-radius ratio and found that in 
our fiducial model, within r < (100 — 200)km the ratio dipped below 0.5, the 
value characteristic of a geometrically thick flow, but only moderately, to 0.3. 
The flow remained relatively geometrically thick at all times. We attributed 
the observed moderate thinning of the accretoin flow to the cooling of the 
flow by the dissociation of helium nuclei into free nucleons. 



limited by the local pressure scale height. 

At radii r ^> nsco, the turbulent dissipation rate (due to 
MRI) is a steeply declining function of radius and this gives 
rise to strong entropy inversion and convective instability. 
The degree of rotational support in the post-shock fluid in- 
creases inward (Lindner et al. 2010). The inner, rotationally- 
supported torus may, according to the Solberg-H0iland crite- 
rion, be convectively stable in the equatorial direction; insta- 
bility is still present in a direction inclined relative to the equa- 
tor, and a fluid element thus transported obliquely, along a 
gyrentrope, eventually mixes with the denser equatorial fluid. 
With this in mind, we develop an effective, spherically aver- 
aged picture in which we postulate that convective heat trans- 
port proceeds according to the prescription of mixing length 
theory (MLT) for a non-rotating atmosphere in which the mix- 
ing length is interpreted as a parameter that hides the com- 
plexity arising from the rotation and vertical stratification. 

The specific angular momentum of the post-shock fluid im- 
mediately following shock formation is only slightly larger 
than that of a circular orbit at risco- The mass of the shocked 
fluid comprises only a small fraction of the mass of the pro- 
genitor star. Barring an extremely steep pre-collapse radial 
gradient in the specific angular momentum of the progeni- 
tor star, the specific angular momentum in the shocked fluid 
varies only over a narrow range of values, and is further ho- 
mogenized by convective mixing. At the smallest radii, vis- 
cous redistribution produces a small, positive gradient in t 
turning over to a small negative gradient at radii at which the 
viscous time (associated with the MRI stress) becomes com- 
parable to the age of the flow. At radii where the viscous time 
is longer than the age of the flow so that no significant viscous 
angular momentum redistribution could have taken place, the 
specific angular momentum is a passive scalar transported by 
convective eddies. We refer the reader to Figure 4 in |Lindner| 
|et al.| (2010), where the near-radial-independence of the spe- 
cific angular momentum of the shocked fluid can clearly be 
seen. 

3. ADVECTION AND CONVECTION 

A fraction of the energy dissipated by the accreting shocked 
stellar envelope is advected into the black hole; the rest is 
transported outward by convection and can power an explo- 
sion. Here we study the competition between advection and 
convection in the innermost segments of the accretion flow 
and attempt to assess the asymptotic luminosity carried by the 
shocked envelope. We work in the spherically averaged pic- 
ture in which fluid variables depend on the spherical radius r 
and represent spheric al av erages over the angular coordinates 
9 and <fi. In Section 3.1 we write relations for the transport 
of inter nal a nd total energy in the inner accretion flow, and in 
Section [3~2l we di scus s the nature of radial force balance in the 
flow. In Section [33] we presen t ou r toy model for the radial 
structure of the flow. In Section [3^4] we justify our adoption of 
a simple equation of sta te that will serve as basis for our toy 
models. In Section 3.5 we review the several key tim e sca les 
characterizing the innermost accretion flow. In Section 3.6 we 



argue that convection is an effective energy transport mech- 
anism only down to some minimum radius; at still smaller 
radii, the energy dissi pate d by accretion is advected into the 
black hole. In Section [3T7| we prov ide e stimates of luminosity 
carried by convection. In Section |3.8|we discuss the impact 
of nuclear processes, and in Section \x9\ we summarize our 
conclusions to help us prepare to undertake an analysis of the 
structure and dynamics of the envelope in Sections [4] and [5] 
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3.1. Energetics and Transport 

At radii much smaller than the radius of the shock, r <C r S h, 
where the accretion flow is in a quasi-steady state character- 
ized by a radial force balance and the bulk inward motion is 
entirely due to viscous accretion, the conservation of internal 
energy can be written in the form 

ds 1 d 2 
v r pT— + -j—[r (F conv + v r pe nuc + F m i x )] = Q visc -Q„, (1) 
dr r L dr 

where v r is the mass-weighted average radial velocity, p is 
the fluid density, s is the specific entropy, e nuc is the specific 
(negative) nuclear binding energy, F com is the heat flux car- 
ried by convection, F m \ x is the flux of nuclear binding energy 
due to convective mixing, Qvisc is the rate of viscous dissi- 
pation that is proportional to the square of the local shearing 
rate, and Q v is the rate of cooling through neutrino emission 
under optically thin conditions. For the purpose of analytic 
transparency and clarity, we have opted not to carry out a for- 
mally self-consistent spherical averaging procedure in which 
all specific-angular-momentum-dependent terms in equation 
(fT) and, depending on the symmetries assumed, the forthcom- 
ing equations could carry additional numerical factors result- 
ing from the ^-dependence of £{r,9, (/>)? 

In the part of the flow where the fluid is in radial force equi- 
librium, the radial motion associated with the viscous angular 
momentum transport occurs with velocity 
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where v is the kinematic shear viscosity and il = l/r 2 is the 
angular velocity. This result applies at the radii that are in vis- 
cous quasi-equilibrium, i.e., where radial angular momentum 
transport rate is approximately independent of radius; note 
also that v r (r) must be continuous and differentiable at any 
local extrema of £(r). 
The conservation of total energy can be expressed as 
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r z pve— + r £ (F com + F mix )}=-Q u , (3) 



where 7 is the adiabatic index of the shocked fluid and $ is 
the gravitational potential. We work in the approximation in 
which the gravitational potential and fluid mechanics are non- 
relativistic; this is clearly not true close to the black hole, but 
for our purposes it will suffice that it be a good approxima- 
tion outside the innermost advective region. In line with our 
hypothesis that the post-shock-formation accretion rate drops 
below the critical value for efficient neutrino cooling, we will 
assume Q v s=s 0. 

In MET, if global compositional gradients and nuclear com- 
position changes inside convective cells can be ignored, the 
heat flux carried by convection is 



1 
Pconv - ~J c pP 
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T ds\ 
Cp dr 1 



3/2 



(4) 



p v^/pj 

3 For example, if the density p is assumed to be spherically symmetric 
and spherical shells are assumed to rotate rigidly, £ <x sin - 8, then the "£" 
appearing in the equations in this section can be interpreted as representing 
2/3 of the maximum, equatorial specific angular momentum, while the terms 
quadratic in £ would require an overall correction factor of 6/5. 



where cp is the specific heat at constant pressure and A con v 
is the convective mixing length. The nuclear binding energy 
is a sum over nuclear species, e nuc = ^■.Ej.JQ/m,-, where £,■, 
Xi, and m,- denote, respectively, the binding energies, mass 
fractions, and nuclear masses of the species. We model the 
mixing of nuclear species in the diffu sion approximation (e.g., 
ICloutman & Eoll|1976l|Kuhfufi|1986| > 



dt 



where 
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(5) 



(6) 



is the velocity of convective cells, and x m i x ~ 1 is a dimen- 
sionless parameter characterizing the efficiency of convective 
mixing. 4 The flux of energy due to the convective mixing can 
then be obtained by multiplying equation d5l> with Ei/rm and 
summing over nuclear species to find 



1 



de„ 



, w u nuc 

^mix — ^rXmix^'conv^convP , ■ 

3 dr 

We assume that the convective motions are subsonic. 



(7) 



3.2. Radial Force Balance and Viscosity 

If the shocked accretion flow has nearly uniform specific 
angular momentum, the fractional contribution of rotation to 
radial force balance decreases radially outward. A relatively 
large mass fraction of the shocked fluid is supported by radial 
pressure gradie nts, and only a very s mall fraction is rotation- 
ally supported (Lindner et al.||20lD) . Let r rot be defined as 
the radius at which the radial pressure gradient force and the 
centrifugal force are equal, -p~ x dP jdr = £ 2 /r 3 ~ jV<i>. We 
will find in the model star that we consider, the black hole 
strongly dominates the gravitational potential at r < 10 4 km 
at all times (see Section E\ and Figure [TJ; at these radii, the 
infalling gas is a negligible perturbation. Thus, here we take 
that V$ ~ GMbhA 2 - Given the weak dependence of £ on 
radius, since the centrifugal acceleration equals a half of the 
gravitational acceleration at r rot , we have 



'"rot 



GM B 



(8) 



The radius r lot should not be confused with the circular test 
particle orbit radius which occurs at £ 2 /(GMbh), again treat- 
ing the potential as Newtonian. While the fluid is in approx- 
imate radial force balance on both sides of r rot as long as the 
latter radius is contained within r s h, the structure of the flow 
changes character at r lot . 

At pressure-supported radii, r S> r Iot , the contribution of ro- 
tation to radial force balance is negligible, and thus a hydro- 
static balance can be achieved if the pressure is P ~ p\ P VQ. 
If the vertical and horizontal pressure scale heights are com- 
parable, i.e., if the flow is thick and radial stratification limits 
the gr owing wavelength of th e MRI, viscosity can be modeled 
with (Thompso n et al.|2005| l 



aXpfl 



(pressure support), 



(9) 



4 For a criticism of the application of MLT to mixing, see e.g., |Ventura et| 
Ial.|dl998l and references therein. 



where Xp = |VlnP| ' is the local pressure scale height. The 
viscous heating rate is Q v i sc = pvo 1 , where a = rdSljdr is the 
shear, With this, the viscous heating rate is 

2 



2vi.sc ~ apX P fl r 



dfl 

~d~r~ 



(pressure support). 



(10) 



Where the flow is predominantly rotationally supported, 
t <C /" ro t, and especially if it collapses into a thin, neutrino- 
cooled disk, in which the vertical pressure scale heigh t limits 
the MRI, we wou ld instead have the thin disk value ( |Shakura| 
' & Sunyaev|l"973) 

v^aP/pQ, (rotational support). (11) 

We will find that the radii of interest are almost certainly unaf- 
fected by neutrino cooling, but the form of viscosity in equa- 
tion ( 1 1 1 is still the appropriate one in the region r < r mt . With 
this, 



Q 



"^'" a n 



dr ) 



(rotational support). 



(12) 



3.3. A Model for Radial Structure 

With the black hole dominating the gravity, V$ sa 
GMsn/r 2 , we model the density, pressure, and specific an- 
gular momentum with power-law profiles, 



p<x.r 



Pocr~ ? , £< 



(13) 



\r, and vertical pres- 
GMftu/r, it follows 



Then, mass continuity d(r 4 v r p)/dr = combined with equa- 
tions (pi and d9l) and hydrostatic balance implies that £ ss 
(5+1 = A + 2. Note that A can be positive or negative. With 
this, for nearly radially-independent specific angular momen- 
tum A ss 0, we have that £ ss 2, so that P ~ ^/oGMbhA and 
v ~ \at. Recall that an approximate radial independence of 
the specific angular momentum is expected given the argu- 
ments we have presented in final paragraph of Section[2]and is 
seen in rotating, two-d imensional numerical simulations (see, 
e.g., |Lindner et al.|2010| Figure 4). 

At rotationally supported radii, r < r rot , if the disk half- 
thickness is one half of the radius, H <~ ■ ' ■ 
sure balance requires (H/r)~ 2 {P/p) 
that P ~ ^pGM Bii /r and from equation (Hi we find that 

v ~ \aGM-RHrJt. Mass continuity and hydrostatic balance 
now imply £ « £+ 1 « 3 — A, Note the slightly different numer- 
ical coefficient multiplying pGM-Qujr in the expression for 
pressure in the rotationally and pressure supported regimes. 

In what follows, we adopt, but cannot rigorously justify, the 
power-law model for the radial dependence of the density and 
pressure in equation ( |13) . The model is certainly ad hoc, but 
it does seem to crudelyapproximate the structure of the solu- 
tions that we have obtained with full, time-dependent, hydro- 
dynamical i ntegrations that we present in a separate, compan- 
ion paper (Lindner et al. |2011[ ), where we find 

£~2, Awl, A«0 (pressure support) (14) 

over a wide range of radii where the shocked fluid is sup- 
ported by pressure. In our analytical toy model, unless we ex- 
plicitly state otherwise, we assume that the power-laws with 
indices in equation ( p"4| i describe the radial structure of the 
pressure supported shocked fluid even where viscous quasi- 
equilibrium has not been reached. The power-law model al- 
lows us to investigate the general properties of collapsar hy- 
drodynamics in the aftermath of the formation of a black hole, 



and illustrate, as we shall see, the potential for an accretion- 
powered explosion, but does not grant us the ability to assess 
the energetics of the explosion with accuracy. 

3.4. The Equation of State in the Innermost Flow 

To allow us to develop an analytical toy model of the in- 
nermost accretion flow, we adopt a simple equation of state. 



Chen & Beloborodov (2007) showed that the the inner, ge- 



ometrically thin, neutrino-cooled disk is on the cusp of de- 
generacy and does not submit itself to r eduction to one of 
the an alytically tractable limits (see, also, Bisnovatyi-Kogan 
|2001 1. However, in our companion numerical investigation 
( |Lindner et al.|20PT| l, we find that following the rapid drop in 
the central accretion rate, the accretion flow no longer cools 
efficiently and is hot and geometrically thick with the densi- 
ties and temperatures reaching ~ 10 8 gcm" 3 and ~2x 10 10 K, 
respectively, at the relevant, innermost radii. Relativistic elec- 
trons and pairs dominate the pressure and the adiabatic index 
remains s=s 4/3. Following a brief transient associated with the 
initial formation of a rotationally supported torus, degener- 
acy remains weak, p, e < kT, within the innermost ~ 1000 km, 
where p, e denotes the electron chemical potential. 

To proceed with our highly simplified analysis of the ampli- 
tudes the various terms in equation <TTJ, we will assume that 
photons and low-density relativistic electrons and pairs dom- 
inate the equation of state in the thick disk so that P oc T 4 
and s oc T 3 /p, and will pretend that the impact of nuclear 
processes on the energet ics c an be ignored (we lift the lat- 
ter restriction in Section |3.8| below). For proton-to-nucleon 
fractions Y e ~ 0.5 and the relevant range of densities and tem- 
peratures, the assumed scalings of the pressure and specific 
entropy are valid as long as (30 /7Tr 6 ){ppc / kT) 6 <C 1, where 
Pp = (3ir 2 p/2m p y> 3 h is the Fermi momentum, or T 3> 1.4 : 



10 10 (p/10 8 gcrrT 3 ) 1/3 K (see, e.g., |Bisnovatyi-Kogan 2001 1. 



In the hot, geometrically thick, low accretion rate regime, this 
condition is satisfied in the inner ~ 1000 km, where we are 
about to argue a critical radius should exist within which con- 
vection will fail to transport the dissipated energy. Our adop- 
tion of a simple equation of state allows us to construct a crude 
analytical toy model for the accreti on flow; for a more a ccu- 
rate numerical treatment, please see Lindner et al. (201 1 1. 



3.5. Time Scales 

For reference, here we provide and comment on several 
time scales characterizing the innermost accretion flow. The 
time scales can be compared with the one on which the energy 
liberated in central accretion accretion flow energizes the stel- 
lar envelope, which in SectionBlwe find is rather long, ~ 10 s 
or longer. We evaluate the characteristic time scales assuming 
the ad hoc power-law radial scalings in equation ( 13 1 with the 
indices given in equation ( [14) , in the regime in which the flow 
is supported by pressure. Iriplace of the dynamical time scale 
we quote the time scale for free fall from rest at infinity 



r 3/2 



TFF' 



VFF (2GMBH)'/ 2 
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3/2 /Af, 



100 km/ 



5 Mr. 



-1/2 



(15) 



The viscous time scale can be estimated by dividing the ra- 
dius with radial velocity, given in equation (p) resulting from 
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angular momentum transport by the viscous torque 



1 



2aAjU 



' 20mS OiT (loknJ 



£ 



10 17 cm 2 s" 



(16) 



Thus, over the relevant ~ 10 s time scale, only portion of 
the shocked, pressure supported material within the inner 
~ 2000km can accrete viscously. Since T conv ex r 5 / 2 while 
r v ; sc oc r 2 with a larger numerical prefactor, we have r conv <C 
r v ; sc at the larger radii in the shocked fluid that convective 
cells have reached. Convection alone will tend to erase gra- 
dients in the specific angular momentum profile, consistent 
with our assumptions that A « 0, and this is clearly seen in 
two-dimensional simulations (see, e.g., Lindner et al 
Figure 4). 



2010 



With the convection we associate the convective eddy radial 
crossing time scale 
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^4ms 



4 r 3/2 

(2GMBH) 1 / 2 
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100 km 
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(17) 



which is only somewhat longer than the free fall time scale as 
a consequence of the steep negative entropy gradient implied 
by our assumed density profile. 

We turn to estimating the time scale for cooling by neutrino 
emission. At densities p ~ 10 7 g cm" 3 (the density dependence 
is relatively weak) and temperatures > 10 10 K, the nuclei are 
almost completely disintegrated and neutrino emission by pair 
capture onto nucleons (Urea) dominates the cooling. At lower 
temperatures, a fraction of the nucleons are in nuclei and neu- 
trino emission by pair annihilation dominates. Ignoring the 
weak dependence of the free nucleon fraction a nd of the Urea 
cooling rat e on density, we c an write (see, e.g., |Popham et al.] 
[1999) |Di Matteo et al.|2002) 



Gcooi~10 25 ergss 'cm 



0.9(r K /10 10 ) 6 (p g cm-Vl0 7 ), 7k> 10 
0.5(7k/10 10 ) 9 , T K <10 



10 



o 'CIS) 



With this, the cooling times are 



7"cool 



(7 _ l)Gcool 

2.2s(r km /lOO) 1 / 2 ( jOgcm - 3 /lO 7 r 3 / 2 (M B H, /5r 1 / 2 

0.6 s (r km / IOO) 5 / 4 (p gcar3 J IO 7 )- 5 / 4 (M BH ,0 /5)" 5 / 4 ; 

where the first version applies at temperatures above ~ 10 10 K, 
the second below. Note for the densities assumed here and 
verified in Section B] and Figure W\ we have that r coo ] 3> 
TVisc, T conv at all radii, which justifies our neglect of cooling 
in the treatment of the flow energetics. We provide the time 
scale for convergence to nuclear statistical equilibrium (NSE) 
in equation ( |53| below. 

An additional time scale is the shock crossing time (see Sec- 
tion[5]below) 



7"shock ' 
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(20) 



where v s h is the shock velocity. Comparing T conv with r s h ck 
for the fiducial choice Mbh ~ 5M Q , we find that for v s h ~ 
10 4 kms _I , the convection can keep up with the shock out to 
~ 600km, while for slower shocks v s h ~ 10 3 kms _1 , convec- 
tion can keep up until the shock has reached ~6x 10 4 km. 
Indeed, we will find that that shock slows down considerably 
after the first ~ 10s. This justifies our assumption that con- 
vection is fully developed in the shocked stellar envelope and 
can transport energy from the innermost accretion flow toward 
the the shock wave in the course of its traversal of the progen- 
itor's interior. The shock, however, must eventually decouple 
and re-accelerate in the outer atmosphere. 

3.6. The Failure of Convection at Small Radii 

We will attempt to compare the relative amplitudes of 
the advection term 2 a( j v = v r pTds/dr, the convection term 

g conv = V • F com = r~ 2 d(r 2 F com )/dr, and the viscous heat- 
ing term gvisc in equation (fill. The detailed form of these 
terms depends on the equation of state; we wor k assum- 
ing that P oc T 4 and s ex T 3 /p (see Section 3.4 1. In the 
pressure-supported regime, with hydrostatic balance imply- 
ing £ m 5+ 1, we find £>adv ex r~ 2 ~^ +x , while the convective 
term is Q com oc r" 3/2 "^, thus the ratio of the advection to the 
convective term is 2adv/2conv oc r" 1/2+A , which for rotation 
laws A < 1/2 increases inward. This opens the possibility 
for the presence of an advection-dominated accretion flow 
(ADAF) at the very smallest radii; at larger radii, we have a 
convection-dominated accretion flow (CDAF; see, e.g.,|Stone 



et al.|1999t|rgumenshchev et al.|2000[rBlandford & Begelman 



2004| l. In the rotationally-supported regime, <2 a dv «r ' 5 A 
and <2adv/<2conv oc r 1/2_A , implying that, to the extent that 
convective transport in the rotationally-supported regime (see 
Section [2]) can be modeled with MET — this assumption can 
only be tested and calibrated with multidimensional hydro- 
dynamic simulations — the smallest radii favor a CDAF. Thus 
for estimating the luminosity carried by the shocked fluid to- 
ward the expanding shock wave, it is key to pin down the 
radius tadaf of the ADAF-CDAF transition and its relation to 
the the radius r rot separating the inner rotationally supported 
region from the outer, pressure supported region. The radii 
~ r rot seem to be the most susceptible to the appearance of an 
ADAF. 

The ADAF regime can occur at radii smaller than a transi- 
tion radius tadaf where the spherically-averaged advection 
fluxes in equation d3), in brackets, exceed the convection 
fluxes F conv and F m i x . This, in general, differs from the cri- 
terion requiring that (Q a dv, Gconv) > Gvisc in the equat orial re- 
gion, which has been employed elsewhere (e.g., Kohri et al. 
|2005| l, since we distinguish spherically-averaged advection 
and spherically-averaged convection. In what follows we will 
conservatively assume that an ADAF is present and 



'"ADAF ^ r rot 



(21) 



so that at r < tadaf, the bulk of the dissipated energy travels 
inward. If ADAF is absent and CDAF extends to the inner- 
most radii, then tadaf in the forthcoming development should 
be replaced by risco; in this case, the convection must carry 
most of the energy dissipated according to equation ( [T2| . 

Comparing the negative (radially inward) radial advection 
luminosity 

2 l p 
^adv = 4-n-r v r r (22) 



7" 



1 



to the positive (radially outward) accretion luminosity 
^acc = 47rr v r p$> 



(23) 

-£+A+l 



in the pressure-supported regime we find that L adv ex r 
while L acc oc r _1 . Again, since £«<5+l = A+2forr> r TOt , we 
find that L a( j v = -2L acc ex r _1 . On the other hand, the convec- 
tion luminosity 



: 4nr z F a 



(24) 



isLo 



-£+3/2 , 



-A-l/2 



In this regime, the energy flux due 



to kinetic energy advection is minus one-half of the viscous 
torque; their sum should be relatively small at the radii where 
L conv dominates energy transport. With A < ~, therefore, there 
is a radius at which the positive convection luminosity L con v 
cannot compete with the inward advection luminosity reduced 
by the outward accretion luminosity L a d v +£ a cc as long as pres- 
sure support holds. Note that the viscous luminosity 



ivisc = -47T7- pvl 



dr 



(25) 



is relatively small compared to L a d v and L acc . 

3.7. Estimates of the Luminosity 

The luminosity that is transported outward from the small- 
est radii through the convective shocked stellar envelope is 
challenging to estimate because of the non-self-similar na- 
ture of the accretion flow. In this section, we consider en- 
ergy transport by hydrodynamic processes: advection, vis- 
cous stress, and convection; we recall that cooling by neu- 
trino emission is inefficient and defer addressing of the role 
of nuclear compositional transformation in the transport of 
energy to the following section. To offset theoretical uncer- 
tainties, we attempt to place multiple constraints on the lu- 
minosity. If the specific angular momentum in the shocked 
region is nearly independent of radius, |A| <§C I, mass conti- 
nuity at the rad ii wh ere dljdt ~ implies 5 f=a 1, As indi- 
cated in Section 3.3 we adopt an ansatz whereby £ = 2, 6 = 1, 
and A~0 throughout the pressure-supported section of the 
shocked region, even the radii where the flow has not reached 
viscous quasi-equilibrium, so that p(r) = p(r ad af)'"adaf / r. 
Since the sum of advection, accretion, and convection lumi- 
nosities, L ac j v +L acc +L conv , must be independent of radius in a 
quasi-steady state, the pressure profile at the radii where con- 
vection does not dominate energy transport should rise inward 
less steeply than radial hydrostatic balance implies; the latter 
is possible, of course, because of the increasing role of rota- 
tion at small radii. 

We hypothesize that L a( i Y (r Iot ) + L acc (r m ) ~ 0, which is in- 
deed satisfied if P(r < r rot ) ~ ^pGM B u/r as we suggested 
above; thus, the net luminosity flowing through this radius, ig- 
noring nuclear and neutrino contributions, is simply the con- 
vection luminosity L = L conv {r rol ). With this, we have that 
L = (L a dv+£ a cc+£conv)('"rot) ~ ^i> where in th e toy model with 
radiation-like equation of state (Section [3 ,4| i, from equations 
(|4) and ([§>, 



vOrot) ~ irGM BH £p(r mt ) 



^convv rot/ 
'"rot 
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(26) 



Note that L\ does not depend on the viscosity parameter a. 
The ADAF-CDAF transition radius tadaf can be defined as 



22k, and d24]i we obtain 



?"ADAF ~ 32a 
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(27) 



which is larger than r rot when A conv (r A DAF) < 2a ' tadaf- The 
advection radius tadaf is larger by a dimensionless factor of 
the form ~ 1 6[A corrv (rADAF)AADAF]~ 2 than the critical radius, 
discussed by |Abramowicz et"aLl ( |2002| , at which the accretion 
velocity |v r | exceeds the velocity of convective cells v conv . 

Combining equation ( 27 1 with equation ( 23 1 we have L = 
(L adv +L acc +L conv )(r A DAF7^ L 2 where 



Lq = L. dcc (r ADAF ) ~ 2iraGM BH £p(r ADAF ). 



(28) 



Since convection can carry the luminosity L to the accretion 
shock, the energy available to power an explosion is also sen- 
sitive to £, A conv , and a; it will be also sensitive to the nuclear 
processes that we have neglected in this toy model thus far. 
Of course, if a is particularly small, neutrino cooling and the 
settling of the flow into a thin disk may start competing with 
advection at the smallest radii. 

An additional lower limit on tadaf can be placed by not- 
ing that this is the smallest radius such that convection can 
transport all the energy dissipated exterior to the radius, 



gvi.sc47rr dr < L com (r) for any r > r ADAF (29) 



which, with equations^, (fT0]>, and ( 24 1 yields the condition 



'"adaf > 2a 2/3 



GM B 



-^conv('"ADAF) 



''ADAF 



-4/3 



(30) 



This differs from our previous estimate only by a numerical 
constant of the order of unity and should be interpreted as only 
a lower limit on the ADAF-CDAF transition because of the 
possibility that advection or nuclear energy transfer through 
disintegration, recombination, and compositional mixing is 
still able to carry inward some of the energy dissipated out- 
side this radius. 

The density at the rotation and the advection radius can be 
related to the immediate post-shock density via / o(r rot )r rot ~ 
P('"adaf)'"adaf ~ p('"sh)'"sh- With this, we have that, for a ~ 
0.1, the two different luminosity estimates are comparable, 

Lz/Li - (8a)" 1 [A conv (r A DAF)/'-ADAF] 2 ~ O(l), which lends 
support to the consistency of our two methods to estimate the 
luminosity carried by the convective envelope. We can thus 
combine equations ( 26 1 and (|28]> to write 



1 f\ 



v conv 

r 



IT (GM BH ) 2 
2 



£ 



r S hp(r S h), 



(31) 



where the dimensionless coefficients in brackets correspond 
to L ~ L\ and L ~ Li, respectively. 

If the gravitational potential is approximately Keplerian, 
hydrostatic balance in the pressure-supported outer parts of 
the post-shock region requires that P ex r~^ with £ = 2, and 
thus, most of the mass and energy resides close to the shock. 
In a medium in which relativistic, low-density electrons and 



the radius at which L a dv =-i C onv From this and equations Q, pairs dominate the equation of state, 4irr 2 F a 



ex r 



-C+3/2 
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hence (8, £) = (1, 2) does not ensure a radius-independent con- 
vective luminosity. In reality, some of the convective lumi- 
nosity is converted into bulk motion with v r > 0, and weakly 
relativistic and degenerate electrons and pairs will dominate 
the equation of state; these effects will reconcile the ra- 
dial momentum equation for the nonrotating outer region, 
v r dv r /dr = — V$- p~ l dP/dr, with the convective flux conser- 
vation d(r 2 F com )/dr = to yield the correct, and likely radi- 
ally dependent, indices S and £. We allow our toy model to de- 
part from local self-consistency; instead, in Section [4] below, 
we require global mass and energy conservation. Beforehand, 
however, we must address the possibility that nuclear compo- 
sitional transformation contributes to radial energy transport. 

3.8. Nuclear Disintegration and Recombination 

Our requirement that neutrino cooling be inefficient at nsco 
can clearly be relaxed to require that it be inefficient at tadaf; 
this is generally the case for relatively large a and central 
accretion rates M <C 0.1 M Q s _1 . Relativistic effects are also 
weak outside r& DAF , especially when the black hole rotates 
rapidly. A more significant concern is the possibility that nu- 
clear composition changes as material passes the shock and 
arrives at tadaf- At densities expected in the vicinity of tadaf, 
which are ~ 10 6 - 10 8 gcm" 3 , and on time scales ~ 0.1 - 10s, 
the principal nucleosynthetic products reach nuclear statistical 
equilibrium (NSE) conditions at temperatures T > 4 x 10 9 K. 
For proton-to-nucleon fractions Y e ~ 0.5, the iron-group-to- 
helium transition takes place at T ~ (5-7) x 10 9 K, and the 
helium-to-free nucleon transition takes place at T > (7- 10) x 
10 9 K. These temperatures can be compared to an estimate 
of the temperature in our model at pressure supported radii, 
r > r roU that reads 



r~1.3xl0 lu KA- a( y/3 

1/4 



GM BH r 
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(32) 



where /3 rad i» is the fraction of pressure in relativistic species 
to total pressure in units of one third. The temperature of 
the accretion flow clearly straddles the temperatures at which 
compositional transitions occur. 

If the helium-to-nucleon transition occurs at radii r,a s com- 
parable to or smaller than tadaf, which can be the case during 
a period following shock formation, then the net energy flux 
through r ADAF should be augmented with 
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where ej^ c w -8MeV/m p is the specific (negative) pre-shock 
nuclear binding energy of the material at the advection radius. 
The second term in brackets in equation p3j ) represents the 
outward energy transport arising from the convective compo- 
sitional mixing. If, with the help of equations Q and |9|, we 
identify 



M = -4-irr pv r ~ 2na£r sh p(r sh ) 



(34) 



with the rate with which shocked material accretes onto the 
black hole, and further identify 



with the rate with which convective cells transport mass radi- 
ally, then we can rewrite equation ( 33 1 as 
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M a 



= 4irr pv c , 



(35) 



It is typical of CDAFs that M conv > M, and thus, depending 
on the relative location of r^ and tadaf, the mixing term may 
reduce the parameter space in which L nuc < and could po- 
tentially even lead to L nuc > 0. 

3.9. Summary on Advection and Convection 

Our goal has been to estimate the rate with which energy is 
carried toward larger radii of the shocked stellar atmosphere 
where the flow has the potential to become "dynamical" and 
engender an explosion. If given such an estimate of the lumi- 
nosity of the inner accretion flow, we are in position to inves- 
tigate the global, time-dependent hydrodynamics of the star 
by treating the central luminosity as a source and by requiring 
global energy conservation. This is the subject of the follow- 
ing section. To attempt to estimate the luminosity, we have 
examined radial energy transport in the inner accretion flow 
which occurs in an approximate quasi-steady-state so that the 
time derivative term in the spherically-averaged energy trans- 
port equation can be neglected and radial force balance is a 
good approximation. Unfortunately, exact analytical treat- 
ment is complicated by the radially increasing dominance of 
pressure support over rotational support in our model. For 
this reason, we have steered away from attempting to iden- 
tify a formal solution to the eigenvalue problem defined by 
the conservation of mass, angular momentum, and energy in 
a flow allowing for advection as well as convection (see, e.g., 
|Abramowicz et al.|2002| |Lu et al.|2 004 who solve the eigen- 
value problem under simplifying assumptions that do not ap- 
ply in the present context). Instead, we have adopted an intu- 
itive approach in which we make educated guesses about the 
radial scaling of the fluid variables. 

Our intuitive approach has allowed us to conclude that 
while the dominant energy transport term at large radii is con- 
vection, at small radii, convection should not be able to com- 
pete with advection, implying that the flow transitions from a 
CDAF to an ADAF at some critical radius. This critical radius 
will generally depend on the viscous stress-to-pressure ratio 
a, the efficiency of the convection (which we parameterize in 
terms of the mixing length), and the angular momentum of 
the accreting fluid, but in our approximate treatment, we are 
not able to pin it down with absolute certainly; instead, we 
provide an approximation in equation ( 27 1 and a lower limit 
in equation OOj, Only the energy dissipated outside the criti- 
cal radius wiirflow outward and contribute to the luminosity 
of the central source. We have obtained two slightly different 
estimates of the luminosity, which we summarize in equation 
pT) , and have also estimated the impact of energy transport 
by nuclear disintegration and recombination in equation p6*| ). 
In the following section, we incorporate the luminosity es- 
timate in dynamical model of the envelope surrounding the 
inner, rotationally-supported accretion flow. 

4. MASS AND ENERGY CONSERVATION 

We proceed to study the global dynamics of the stellar en- 
velope outside of the innermost, rotationally-supported accre- 
tion flow. A fraction of the envelope has passed the accre- 
tion shock; the kinematics of the shock, as we shall see, is 
the principal determinant of the conditions in the innermost 



flow, while the luminosity transported outward from the in- 
nermost flow is the driver of the shock's dynamics. Keeping 
in mind this interdependence, we model the density and radial 
velocity structure of the fluid flow inside and outside the shock 
radius and require mass and energy conservation. Normally, 
mass, momentum, and energy conservation are imposed in the 
form of jump conditions applied to the shock transition itself, 
but this approach generally leads to violation of global mass 
and energy conservation, unless one solves the full, time- and 
radius-dependent hydrodynamic transport equations. For the 
purpose of analytic transparency, we opt for a simpler ap- 
proach in which we adopt a model density and velocity field, 
and then require global mass and energy conservation. This 
may violate local conservation at the shock. 

If at mass coordinates just exterior to those first passing 
through a nascent accretion shock, the stellar pre-supernova 
density profile is approximately p*,,=o ~ M+jA-Kr+r 2 , where 
M+ the stellar mass enclosed within some radius r*, then 
the subsequent flow of material into the shock may resem- 
ble the self-similar collapse of nearly-hydrostatic isothermal 
spheres of Shu (1977i. This solution is characterized by a 
critical point (a rarefaction front) that recedes into the pre- 
collapse envelope at the speed of sound. The critical point 
starts traversing the star outward at t = 0, and thus it is well 
ahead of the shock front, which starts traveling after tens of 
seconds, even if the shock travels supersonically. The density 
profile of the collapsing envelope in the self-similar solution 
steepens from p* ex r" 1 in the immediate vicinity of the crit- 
ical point to p+ ex r" 3 / 2 at radii r <C r cn \ ~ (5 — 10) X 10 4 km, 
where the critical point traverses the star at the sound speed 
dr cnt /dt = (7*GM^/2r*) 1 / 2 , where 7* is the adiabatic index 
of the stellar envelope. We adopt p* ex r" 3 / 2 for the in- 
falling envelope. The fluid velocity can be approximated via 
v* ~ (1 - r/ r crit)vff> where Vff(r, t ) = -[GM(r, t)/r] '/ 2 is the free 
fall velocity. 

Let M t be the stellar mass enclosed within some radius r*. 
Then we model the stellar density profile via 



' M sh /(2^r 2 h r), 



r < r sh , 



P(r) = { Mj{4KrS^I 2 ), r sh < r < r cnt , 



.M*/(47nv- 2 ), 



<r<r±, 



where 



M sh = 27rr sh p(r sh ) 



(37) 



(38) 



is the mass of the shocked fluid. The equation of mass conti- 
nuity in the post-shock fluid, dp/dt + r~ 2 d(r 2 v r p)/dr = 0, im- 
plies, for the outer region r ~ r s h in which viscous accretion 
can be ignored 

d . r sh 



dt 



■In 



M. 



1/2' 



(39) 



and with this we model the radial velocity profile via 



f rd\n(r sh /MH 2 )/dt, r < r sh , 

v r {r)= I -(l- r / rcrit )[GM(r,0/r] 1 / 2 , r sh < r < r cr ,. (40) 

u 



'cnt> 

rent <r<r+, 



where M(r, t) = M B h + J 4-irpr 2 dr. Finally, in the outer parts 
of the shocked region, r ~ r s h, the momentum of the shocked 
fluid contributes negligibly to radial pressure balance, and 
thus the shocked region is quasi-hydrostatic. We model the 



pressure profile via 



r < r sb , 



P(r) = < \ [prV$];? (r/r cnt )"T* , r sh < r < r crit , (41 ) 



i[prVfc]' 



r c m <r<r+, 



where, again, $ is the total gravitational of the black hole and 
the mass distribution in equation ( 37 1, 
Mass conservation requires that 



dt 



M B H+Af sh + / \-Krpdr) =0 



'Mi 



(42) 



The mass of the black hole at the moment of shock formation 

is 

1 »*crit(*sh) , 



3 r* 



-M„ 



(43) 



where Mbh.o is the initial black hole mass immediately fol- 
lowing core collapse in excess of the mass available in the 
core from the inward-extrapolated profile p+ ex r~ 2 . Since the 
accretion rate into the black hole just prior to shock forma- 
tion M(t < f s h) ~ ^(M+/r+)dr C nt/dt is similar to the accretion 
rate into the shock after shock formation, and the black hole 
accretion rate experiences a drop as the shock starts traveling 
outward, we expect dM^jdt ^> dMsu/dt. 

Pretending that r* ~ 10 5 km is the true edge of the star and 
that efficient neutrino cooling, if anywhere, occurs only at r < 
r\DAF, global energy conservation can be written as 



L+L auc = -(U + K+W + E nuc ) 
dt 



d 
dt 



GMv, 



- I Airpr dr, 



(44) 



where L is the net energy flow rate at tadafj while U, K, W, 
and £n UC , are the total internal, kinetic, gravitational poten- 
tial, and nuclear energies in the annulus tadaf < r < r+. The 
contribution of nuclear composition change to the energy flux 
through the advection radius, L nuc , is significant only when the 
helium-to-free nucleon disintegration occurs at radii > tadaf- 
After the shock has expanded far beyond the disintegration 
radius, we can assume that the mass in free nucleons evolves 
very slowly in time, and thus, dE nuc /dt f=a 0. 

The total internal energy is the sum over the shocked and 
unshocked regions 



U 



7-1 



-4irrdr+ 



7* 



■1 



-47rr dr, 



(45) 



where the pressure is calculated following the model in equa- 
tion (41 1. The kinetic energy is calculated via 



K = 



1 



-pv r 4Trr dr. 



(46) 



The total gravitational energy W of the density in equation 
( [37] i in the presence of self gravity and the gravity of the black 
hole is straightforward to calculate, 



W = 



-9^ + U p )p4nr 2 dr, (47) 

r 2 



where the gravitational potential of the fluid outside of the 
black hole is related to the density via 4irGp = V 2 $ p . 



10 



SUPERNOVAE POWERED BY COLLAPSAR ACCRETION 




slowly in time, we do not carry out the formal integration and 
instead approximate 



(cm) 



FIG. 1 . — The densit y p, r adial veloc ity v r , and pressure P of the toy model 



as given by equations 1371, (401, and 141), for M k = 12. 5Mq, r+ 



10 iu cm, 
7* = 1 .4, Mbh.o = 2.5M Q , £ = 10" cm 2 s" 1 , a = 0. 1, A conv /r = 0.75, L„ uc = 0, 
( s h = 20s, and t—t S h = (0 .5, 1,2, 4, 8) s; the convective luminosity was calcu- 
lated from equation (28). The radial velocity does not take into account the 
negative radial velocity resulting from viscous angular momentum transport. 
The density jump at the shock is smaller than it should be because our adopted 
density profile for the infalling envelope, p <x r~ 3 ' 2 for r s h < r < r cr ; t is too 
steep at radii just smaller than r cr ; t , where p <x r _1 . 



5. SHOCK DYNAMICS 



Equations (42 » and (041 can be construed as constraints 
relating the shock velocity to the shock radius and time, 
Vsh(/sh,0- The equation dr^/dt = v s h can then be integrated to 
solve for r s h(f )■ Since we expect that the shock velocity varies 



r s h~(f-fsh)vsh, 



(48) 



where, as before, ? s h denotes the shock formation time. This 
allows us to estimate v$h(t) and r s h(f)- In FigurefT] we show a 
typical time evolution of the density, radial velocity, and pres- 
sure profile for a fiducial stellar model with M+ = 12.5M , 
r*= 10 10 cm, 7*= 1.4, M BH ,o = 2.5M Q , and I = lO^cmV 1 . 
This model approximates the density profile of the fully 
mixed pre-supernova Wolf-Rayet model 16TI of Woosley & 
Heger (2006) in the radial range 2 x 10 8 cm < r < 5 x 10 y cm. 
At radii r > jr+, our toy model overestimates the stellar enve- 
lope density, which declines increasingly steeply with radius, 
and thus it underestimates the shock velocity after the first 
~ 10 s. Note that the black hole dominates the mass enclosed 
f or r -C 4 x 10 4 km at all times, which justifies our neglect of 
the infa lling mass in the calculation of the force balance in 
SectionEHl 

The density jump at the shock in Figure [TJ is substantially 
smaller than it should be; e.g., in the strong-shock limit — 
which is not always reached here — the density jumps seven- 
fold for 7=1. The anomaly seems to be an artifact of our as- 
sumption that the density in the region r s h < r < r cr j t is a pure 
power law p oc r~ 3 ' 2 . In reality, as illustrated in the Shu Q1977 i 



solution for the self-similar collapse of isothermal spheres, the 
density profile is less steep, p oc r , just inside r cr j t , and thus 
our immediate pre-shock density is an overestimate. Since 
the immediate post-shock density is calculated independently 
from global mass conservation, the ratio of the two densities 
as seen in Figure[TJis also an overestimate. In spite of this, the 
model does conserve mass and energy globally. 

For the viscous stress-to-pressure ratio, we adopt a = 0.1 
and in Figure |2l we show the evolution of the shock velocity, 
the total energy of the stellar envelope, and the rate with which 
the envelope is accreting onto the black hole for three val- 
ues of the convective mixing length, A conv /r = (0.5,0.75, 1.0), 
with and without nuclear disintegration losses. The lumi- 
nosity carried by the post-shock region was calculated using 
equation ( 28 • so that L = L-i (the second case in brackets in 
equation 3T ) and thus, v s h(f) and E(t) depend only on the 
ratio (A CO nvA) 4 /(a^)- 

From its pre-shock value of M(t < t sh ) m O.14M s _1 , the 
accretion rate has dropped to M = (0.003 -0.03)M Q s" 1 after 
the first second from shock formation, and to M = (0.001 - 
0.004)M Q s _I after ten seconds. The steeper drops occur 
in the more energetic shocks with larger convective mixing 
lengths; at such low accretion rates and a ~ 0.1, neutrino 
cooling is negligible compared to viscous heating, at least at 
r > /"adaf ^ 0.5 x 10 7 cm; this justifies our leaving out of the 
cooling term in equation ( |44| >. If the gamma ray luminosity 
of the LGRB prompt emission is controlled by the rate with 
which material is accreting onto the black hole, then the steep 
drop in accretion rate associated with shock expansion could 



explain the termination of the prompt emission ( Lindner et al. 
12010) . 

The model with A con v/'" =1.0 acquires positive energy, and 
the envelope is unbound and poised to explode at ~ 5 s after 
shock formation, when the shock velocity is ~ 10, 000km s~ . 
The model with A conv /r = 0.75 is on track to acquire positive 
energy after the shock reaches r cr i t ~ r*, which are shock radii 
that our toy model is not designed to handle. The model with 
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Aconv/r = 0.5 does not appear to evolve toward a globally un- 
bound state, though of course, for a realistic pre-supernova 
density profile that declines steeply with radius at r > ^r+, 
the shock will ultimately emerge from the star and unbind a 
fraction of its mass. The latter model is sensitive to nuclear 
disintegration losses; with L nuc sw -8MeVM/m ;) , the shock 
seems to stall at ~ 2 x 10 9 cm. 

These results indicate that the potential for explosion in 
collapsar-accretion-powered objects depends critically on the 
efficiency of convection. We are not aware of a numerical cal- 
ibration of the effective convective mixing length A conv , if the 
latter is defined as the mixing length that gives a MLT heat 
flux equal to the true energy flux carried by convection in the 
regime, characteristic of supernovae, in which the convective 
velocities are comparable to the sound speed. Such a calibra- 
tion would improve the toy model presented here. 



The 2.5D axisymmetric hydrodynamic simulations of Lind 



ner et al. ( 2010|i, which were carried out w ith a realistic equa- 
tion of state (Tirnmes & Swesty 2000|l, developed a fully 



convective flow in the shock downstream. These runs were 
restricted to the domain with cylindrical radii R > R m [ n = 
(0.5-2) x 10 3 km and the stress-to-pressure ratio was a « 
0.01. The luminosity carried by the convective envelope was 
^conv ~ 0.05 c s P, where c s is the sound speed, and the shock 
velocity is Vgh ~ (0.5- 1.5) x 10 3 kms . If we artificially set 
?"adaf = ^min an d chose a = 0.01 and f~3x 10 17 cm 2 s" 1 , then 
o ur toy model reprod uces the relatively low shock velocities 
(2010). An additional complication not in- 



Lindner et al. 



in 

vestigated in Lindner et al. is the exothermic and endothermic 

compositional change in rising and sinking convective cells. 

6. NUCLEOSYNTHESIS AND NICKEL 

The production of an optically bright supernova requires 
the synthesis of a substantial mass of 56 Ni. This requires that 
a substantial mass of the shocked stellar envelope be heated 
to temperatures >5x 10 9 K. Also, the reprocessed mate- 
rial must freeze out into iron group elements. Finally, the 
proton-to-nucleon ratio during freezeout must be Y e s=s 0.5. To 
check whether nickel may indeed be synthesized in accretion- 
powered explosions, we will examine these requirements, re- 



spectively, in Sections 6.1 6.2 and 6.3 but first, we briefly 



review some of the different scenarios 

In the standard model for core collapse supernovae, nickel 
is synthesized when the shock is fast and the immediate 



post-shock temperat ure is suf ficiently high (e.g., Woosley 
& Wea ver 1995; He ger et al~||2003| ). Nucleosynthesis cal- 



culations for LGRB supernovae and their ultra-energetic 
the hypernovae — typically emplo y a piston to ac- 
celerate supernova ejecta (e.g., [Maeda &~T ominaga |2009[ 



Dessart et al. |2011| l, or inject an energetic jet (e.g., |Tomi-| 
naga et al. 2007), or apply heating in the downstream of the 



stalled shock (e.g., |Fryer et"aL 2006) to initiate an explosion. 
These studies find that high nickel masses inferred in the su- 
pernovae associat ed with LGRBs require the injection of en- 



ergies > 10 52 erg ( Tominaga et al.|2007 Maeda & Tominaga 



2009). However, the physical mechanism that deposits such 
large energies in the stellar envelope remains to be elucidated. 
In a different scenario, nucleosynthesis in the collapsar sce- 
nario occurs in an freely expanding outflow, fireball or wind, 
coming from a disk of material accreting onto the black hole. 
For the outflow to synthesize 56 Ni, the inner disk must be 
nondegenerate so that proton-neutron equality can be main- 
tained, which is possible with moderate accretion rates; al- 
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FIG. 2. — The shock velocity v 8 h, the total energy of the fluid outside the 
black hole E, and the accretion rate onto the black hole M = dM%\\jdt as 
a function of time from shock formation, t-t s h, for the fiducial model with 
M„ = 12.5M©, n = 10 10 cm, 7* = 1.4, M Bli , = 2.5M Q , £ = 10 17 cm 2 s _1 , 
a = 0.1, and fa = 20 s; the convective luminosity was calculated from equa- 
tion (28) . The solid curves are ignoring nuclear disintegration, L nuc = 0, 
for Acorn/i" = (0.5,0.75,1.0); larger convective mixing length give faster 
shocks, more energy deposition, and lower accretion rates. The dotted lines 
are the same but with maximum energy loss due to nuclear disintegration, 
inuc = -8MeVM/m ;) . The pre-shock accretion rate, not shown in the figure, 
isM(f <f sh )«O.14M s- 1 . 
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ternatively, neutron-proton equality can be re-established in 
the wind, and simultaneously, the material must not fr eeze 
out in the expanding wind too quickly to produce nickel (Be- 
|loborodov|2 003; Pr uetet al.|20 03, 2004). These possibil ties 
are clearly very interesting, but they require the presence of 
an open funnel through which the outflow from the inner ac- 
cretion disk can escape. 

A funnel-like density distribution is undoubtedly present at 
small radii where rotational support is competitive with pres- 
sure support. It is not clear, however, that the funnel can be 
open at somewhat larger radii, where rotational support is not 
significant. The axial region may further be overpressured by 
the hot cocoon produced during the electromagnetic jet's first 
traversal of the star. Lacking a funnel, the disk outflow en- 
counters infalling stellar layers. In this regime, however, the 
mechanics of nucleosynthesis in the collapsar must be exam- 
ined in the context of the interaction and mixing of the outflow 
with the (shocked) stellar envelope. To attempt to understand 
the implications of the interaction of convection-like outflows 
from the hot inner region with the cooler, but more massive 
layers of the shocked stellar envelo pe, we ad apt some of the 
useful approximations developed by Beloborodov (2003) and 
|Pruetetal.| ( T2"004l i. 

6.1. The Mass Reprocessed to NSE 

While only a small fraction of the shocked fluid is hot 
enough to allow nuclear burning, because of the pervasive 
convection in the shock downstream, a much larger fraction 
of the shocked fluid can be circulated through the hot inner 
region and can thus be reprocessed into free nucleons, he- 
lium, or the iron group elements. Let tnse denote the radius 
within which NSE among the principal nucleosynthetic prod- 
ucts is established on a convective eddy crossing time. We 
would like to calculate the mass fraction of the shocked fluid 
that under the action of the convective mixing visits the radii 
t < '"nse- For this, we must solve equation 0\ with X denot- 
ing the mass fraction of unreprocessed elements subject to the 
boundary condition that none of the unreprocessed elements 
survive inside r N $ E . We seek a quasi-steady-state solution 



d_ 

Or 



47rr 



1 



j A,. 



dX 



n, Xmix^conv^convP ^ 

3 or 



= 0. 



(49) 



the shocked fluid that is reprocessed through tnse then equals 



fi 



M x (r NS E) r sb 



NSE' 



M s 



sh 



Vsh 



I / ^NSE 

3 V r sh 



1/2 



AconvvshJ ^convv sh/ 



''sh 



»'sh 



(52) 



Khokhlov (1991) approximated the time scale for conver- 
gence to NSE via 



T NSE ' 



^c m -3exp[179.7/(7k/10 9 )- 



■40.5] s, (53) 



which is co nsistent with the more recent estimate of ICalderl 
|et al.| ( |2007) . Setting t NS e ~ 0.1 S, this yields r N sE ~ (0.5- 
2) x 10 s cm where the temperatures are 7nse ~4x 10 9 K. The 
typical total reprocessed mass prior to the final acceleration of 
the shock (for r s j, < rj for our fiducial toy model is 



M NSE -/Ns E M sh . 



(54) 



When at time U ~ '*/v s h('+) the shock radius reaches the edge 
of the stellar model, r s h ~ r+, the reprocessed mass becomes 



M^ s |~O.15M 0Xl 




V10M, 



(55) 



In deriving equation p5| ), we have assumed that the convec- 
tive velocity is of the form given in equation (|50j>. 

It is possible that substantial additional nuclear reprocess- 
ing inside r N sE can take place over a longer period after the 
shock has proceeded to accelerate down the steep density gra- 
dient of the outer stellar envelope and break out of the star. 
Thus, the NSE mass estimate quoted in equation ( 55 1 can be 
considered a lower limit. If following shock breakout at ~ f* 
the density inside the original stellar radius r* decreases ex- 
ponentially, e.g., on a free fall time p t>h ex exp[-Vff(r + )f /r star ], 
then following shock breakout an additional mass is repro- 
cessed through NSE and can be estimated via 



In MET, the convective velocity can be estimated from equa- 
tion Q, which, assuming the power-law pressure and density 
profiles, p ex r~ 6 and P ex r~* and the simple equation of state 
discussed in Sectionl3~4lreduces to 



V ™~\\\t-\ S 



1/2 



A r 



v ff! 



(50) 



where vg is the free fall velocity. At small radii, for the 
purpose of a rough estimate, we can assume that vg ~ 
(GMbh/V) 1//2 so that v conv ex r~ x l 2 . If, as before, 5=1 and 



Aconv oc r, we can rewrite equation ( 49 1 as 



K^f)-- 



(51) 



™NSE 



M x (t)dt 



r* 



Vff(r*) 

Vsh(r*) 
Vff(f±) 



M x (t+)exp 



Vff(r+)t± 



exp 



Vff(r+) 
v S h(r*) 



M' <h 



(56) 



This shows that if the average shock velocity inside the star 
is large compared to the free fall velocity at the stellar sur- 
face, v s h 3> Vft, most of the reprocessing to NSE takes place 
after the shock leaves the star. To arrive at this conclusion, we 
have employed a number of extremely crude approximations; 
a more accurate approach would clearly require carrying out 
a time integration on a realistic model star. 



where = ^-k\ 



independent coefficient. 
X(r NSE ) = 0mdX(r sh ) = 



mi\ ' , >nvA C onv'' 1 ' 2 p is an approximately radius- 



Integrating this twice and setting 
1, to obtain the mass flux of species 



X through r N sE is M x (r NS E) = 2 0r NSE- The fraction / NSE of 



6.2. Freezeout into the Iron Group 

What fraction of the reprocessed fluid can turn into iron 
group elements? Investigations of the nucleosynthetic foot- 
print of freely expanding winds have been carried out by many 
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authors (e.g., Beloborodo v|2003j [Pruet et al.|2003||2p04| [Na 

gataki et al. 2006; Surman et al. 2006; Maeda & Tominaga 



2009, Metzger 2011), but we are not aware of a systematic 
investigation or nucleosythesis in a quasi-hydrostatic, con- 
vective atmosphere stradling a region in NSE and a frozen- 
out atmosphere. We anticipate carrying out multidimensional 
simulations of turbulent convection in the presence of nuclear 
burning to learn about the compositional yields in such flows. 
Here, we attempt to harness the expanding wind solutions 
by applying them to individual convective cells. This is un- 
doubtely extremely crude, but is consistent in spirit with the 
n ature of the a pproximations entering the derivation of MLT. 
|Pruet et aL] ( 200 4 ) calculated the mass fraction of the iron 
group (more precisely, of 56 Ni since they assume mildly 
proton-rich conditions, Y e = 0.51) in a freely expanding col- 
lapsar wind as a function of the entropy per baryon in 
units of the Boltzmann constant, S = (m p /k B )s and the 

variable M w j n d/vJ, ind . Defining the dimensionless parame- 
ter n = (M wind /v^ nd )/[O.lM s-7(O.lc) 3 ], an X Fe > 50% 
freezeout into the iron group requires /i > (2, 13,40) for S = 
(20,30,40), while, similarly, anX Fe > 25% freezeout requires 
(j, > (0.1,0.6,2,4) for S= (20,30,40,50). The freezeout into 
iron group prefers low entropies, high densities, and slow con- 
vection. 

We can apply the [Pruet et al.| (2004) result to a single rising 
convective cell by identifying M wind /v^ ind with 47rr 2 / o/v^ onv . 
In the absence of degeneracy, the entropy per baryon at the 
helium-to-iron group boundary can be estimated from equa- 
tion (12) in Pru et et al.| ( |2004) , 



T 3 1 

S^5.21 MeV +- 

■-• 4 



5 /10 8 



,3/2 



15.4+ln 



A> 



,/10 8 



(57) 



which, with T = T NSE - 4 x 10 9 K, becomes S~ 21/p 6 +4.6- 
0.251np6, where p(, = p(r N sE)/10 6 gcm~ 3 . Electron de gen- 



eracy sets in at temperatures r deg <2x 10 9 p6K (e.g., Be- 
loborodov 2003), below the temperature at which the iron 
group freezeout occurs. If the velocity of convective cells is as 
given by equation ( 50 1 with the exponents given in equation 
( |14) , we obtain 

3 r x / x -\ -2 

. r %Pb 



/'' 



Ms 



A C onv( r NSE) 

'"NSE 



(58) 



where r& = rNSE/10 8 cm, and Ms = Mbh/5M q . This esti- 
mate, which may be excessively conservative, suggests that 
efficient freezeout into iron group elements requires /cK/nse) ^ 
10 6 gcm" 3 . Such densities are clearly realized in the convec- 
tive accretion flow, but are probably not realized in the rela- 
tivistic (or nearly relativistic) axial jet where freezeout is into 
a-particles. 

6.3. Nickel Synthesis and Implications for Supernovae 

The proton-to-nucleon ratio of the stellar envelope entering 
the accretion shock is Y e « 0.5, and with this value, 56 Ni dom- 
inates the composition of the iron group products produced in 
the convective accretion flow. However if significant delep- 
tonization operates at ~ few x risco, convection may transport 
the neutron-rich fluid near ~ r N $E and thus tip the balance in 
favor o f iron and the lighter iron group isotopes. Beloborodov 
(2003) derives an estimate of the equilibrium value of Y e (T,p) 
in at most mildly degenerate matter that is transparent to neu- 
trinos and applies it to a rotationally-supported accretion flow 



with vertical scale height H ~ |r to conclude that Y e drops 
below proton-neutron equality when accretion rates exceed 
M >M„ where 



M„ = O.O55M s" 



(£) 



r„ 



1/2 



M 



BH 



5 Mr. 



(59) 



where r g = 2GMbh/c . Our toy model suggests that the ac- 
cretion rate drops well below M„ very quickly following the 
initial shock formation, and this implies Y e > 0.5, where, at 
densities p ~ 10 8 gcm~ 3 characteristic of the innermost disk, 
56 Ni dominates the iron group. 

These crude estimates make us optimistic that supernovae 
powered by collapsar accretion can synthesize nickel masses 
similar to those required to explain the light curves of su- 
pernovae associated with LGRBs. More detailed work is re- 
quired to characterize the interplay of convection and nucle- 
osyntesis in the shocked, pressure-supported accretion flow 
of a collapsar. A prediction of the present model is that in 
the supernova ejecta, nickel is mixed with a much larger mass 
of unburned stellar material. This mixing produces a super- 
nova with a steeper initial rise that is brighter at early times 
than a spher ically-symmetric explosion (see, e.g., |Woosley & 
Bloom 2006). While our toy model assumes a quasispherical 
Shockwave, the global structure of the explosion should be- 
come aspherical just prior to and following shock breakout, 
and certainly on time scales of ~ 100 s, with higher entropy 
material outflowing near the rotati on axis, as seen in the ide- 
alized 2.5D simulations of [Lindner et al.| ( |2010| l. Su ch as- 
phericities can be detected through spectr opol arimetry (Wang 
& Wheeler 2008) and spectroscopy (e.g., Tanaka et al. 2009) 
and if present at shock br eakout can also be inferred from the 
breakout light curve ( |Couch et al.|20 1 1 ). We will address the 
structure of the ejecta and the implications for the supernova 
light curve and other observational properties elsewhere. 

7. CONCLUSIONS AND DISCUSSION 

With the aim of shedding light on the mechanism that pro- 
duces Type Ic supernovae in LGRB sources, we have devel- 
oped a toy model for the accretion of a rotating stellar enve- 
lope onto a black hole in the aftermath of stellar core collapse. 
The purpose of the toy model is to test the ability of collapsar 
accretion to produce supernovae, and identify aspects of the 
problem, such as the nature of the ADAF to CDAF transition 
and the mechanics of convection, that require further investi- 
gation. 

The spherically-averaged toy model for a rotating collaps- 
ing star assumes that no prompt explosion prior to black hole 
formation takes place. The model is constructed to globally 
conserve mass and energy. We track the dynamics of the out- 
ward traveling shock wave that forms when infalling stellar 
layers have sufficient angular momentum to be held up the 
centrifugal barrier and circularize around the black hole. 

The shocked fluid, heated by the dissipation of MHD tur- 
bulence that we model with a viscous shear stress term, is 
convective; we treat this convection in the mixing length ap- 
proximation. Some of the dissipated energy is advected into 
the black hole; the rest is transported by convection through 
the expanding shocked region and is available to power a su- 
pernova. The amount of energy delivered to the stellar enve- 
lope depends on the location of the boundary of the ADAF at 
small radii and CDAF at large radii. 

The ADAF to CDAF transition is particularly sensitive to 
the effective convective mixing length. If the mixing length 
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is sufficiently large, our model can acquire positive total en- 
ergies of at least ~ 10 51 erg over the course of 10 seconds 
or longer, which lays open the possibility of a supernova. It 
does not seem, however, that the mechanism could produce 
"hypernova"-like energies (> 10 52 erg), at least not in Wolf- 
Ray et progenitors. 

The rate with which shocked stellar fluid accretes onto the 
black hole drops drastically following the inception of the ac- 
cretion shock, and thus, losses to neutrino emission are negli- 
gible. If the luminosity of the LGRB prompt emission is cor- 
related with the accretion rate, then the abrupt termination of 
the prompt emission and the steep decline of the early X-ray 
afterglow can be interpreted as a consequence of the accretion 
rate drop (see, also, Lindner et al. 2010). 

Because of the rapid convective mixing, tens of percent of 
a solar mass can be reprocessed through the hot inner radii 
of the accretion flow where NSE is reached on a dynamical 
time. Some reprocessing takes place after the accretion shock 
breaks out of the star. Conditions are favorable for the freez- 
ing out of the reprocessed matter into 56 Ni, which due to pre- 
shock-breakout convection should be intermixed with a much 
larger mass of hydrostatic a-elements in the stellar ejecta. 

We have assumed throughout that the specific angular mo- 
mentum of the initial star increases more or less mono- 
tonically outward, as one mig ht expect in fully-mixed pre- 
supernova models (Woosley & Heger 2006). If this is not 
the case, then the nonmonotonicity (see, e.g., |Heger et al. 



of a stellar core into a black hole occurs. It would be interest- 
ing to check whether the collapse in to a black hole in the core 
of a rotating supermassive star (e.g., Fowler 1966; Bond et al. 
1984; Fuller et al. 1986; Baumgarte & Shapiro 1999; Shibata 
& Shapiro 2002), could, as we find here for WR stars, lead to 
an unbinding of a significant fraction of the remaining stellar 
envelope. 
In an attempt to elucidate the rapid formation of massive 



black holes in early galaxies, Begelman et al. (2006, 2008) 
and |Begelman ( 2010|> have proposed that a black hole can 



2000 2005 Petrovic et al. 2005 ) might have interesting con- 
sequences tor the evolution of the accretion rate. For exam- 
ple, the accretion rate may surge if the average specific angu- 
lar momentum in the shocked region drops below the critical 
value for rotational support near ISCO, and this might result 
in a "flaring" in th e LGR B X-ray light curve (se e, e.g., [L6pez^ | 
Camara et al .2010; Perna & MacFady 



Wolf-Rayet stars that seem to be the most plausible LGRB 
progenitor candidates were the target this inquiry, but the 
analysis can be adapted to other contexts in which the collapse 



form at the center of a large accumulation of gas (~ 10 M Q ) 
in a gas-rich primordial galaxy. The black hole subse- 
quently accretes the gas at the center of the rotating, pressure- 
supported gaseous object, a "quasistar," in such a way that the 
object remains gravitationally bound. In this picture, the ac- 
cretion rate settles in a quasi-steady state in which the energy 
dissipated at the innermost radii is transported convectively to 
the outer radiative zone; the latter, thanks to an internal self- 
regulating adjustment in the structure of the quasistar, carries 
a radiative flux that remains below the Eddington limit. Our 
results suggest that the fate of the quasistar may depend on 
the energetics of the relatively short period in the immediate 
aftermath of black hole formation and that the path to a self- 
regulating quasi-steady state deserves further inquiry. 
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